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Abstract. The QCD critical point can be found in heavy ion collision experiments via the 
non-monotonic behavior of many fluctuation observables as a function of the collision energy. 
The event-by-event fluctuations of various particle multiplicities are enhanced in those collisions 
that freeze out near the critical point. Higher, non-Gaussian, moments of the event-by-event 
distributions of such observables are particularly sensitive to critical fluctuations, since their 
magnitude depends on the critical correlation length to a high power. We present quantitative 
estimates of the contribution of critical fluctuations to the third and fourth moments of the 
pion and proton, as well as estimates of various measures of pion-proton correlations, all as 
a function of the same five non-universal parameters. We show how to use nontrivial but 
parameter independent ratios among these more than a dozen fluctuation observables to discover 
the critical point. We also construct ratios that, if the critical point is found, can be used to 
overconstrain the values of the non-universal parameters. 



1. Introduction 

One of the distinctive features of the phase diagram is the critical point at which the first-order 
transition between hadron matter and QGP ends. We currently do not have a systematic way 
of locating this point from first principles as model and lattice calculations face many challenges 
and much work still needs to be done in order to overcome them. (For review, see e.g., pQ.) 
In the meantime, if the critical point is located in a region accessible to heavy- ion collision 
experiments, it can be discovered experimentally in current and future experiments at RHIC, 
BNL and SPS. It is therefore important to define and select experimental observables that will 
allow us to locate the critical point, if it is located in an experimentally accessible region. 

Upon scanning in the center of mass energy yfs and thus in fig, one should then be able 
to locate (or rule out the presence of) the critical point by using observables that are sensitive 
to the proximity of the freeze-out point to the critical point [21 [3], such as pion and proton 
multiplicity fluctuations. As we vary y/s, if the freeze-out point approaches the critical point, 
we would see an increase and then a decrease, as we move away from it, in the fluctuations in the 
number of these particles. In this paper we describe how to use this non-monotonic behavior in 
the fluctuations of particle numbers near the critical point as a probe to determine its location. 



1 Talk given by C. Athanasiou at Hot Quarks 2010. 



2. Results 

2.1. Cumulant calculations 

The critical mode, a, is the mode which develops large long wavelength fluctuations at the 
critical point. Its screening mass is related to the correlation length £ by m a = In the case 
of an infinite system, £ diverges at the critical point. In reality, £ reaches a maximum value at 
the critical point but does not diverge because as it cools the system spends only a finite time in 
the vicinity of the critical point. Estimates of the rate of growth of £ as the collision cools past 
the critical point suggest that the maximal value of £ that can be reached is around 1.5 — 3 fm 
[UEIE], compared to the natural ~ 0.5 fm away from the critical point. Near the critical point, 
the a 3 and <r 4 interaction couplings are given by A3 = A3T (T£) _3//2 , and A4 = A4 (T£) _1 , where 
the dimensionless couplings A3 and A4 are universal and they have been determined for the Ising 
universality class [7]. Throughout this paper we shall use A3 = 4 and A4 = 12 as benchmark 
values, because these are the midpoints of the ranges of values known for these constants IB]- 

The long wavelength fluctuations in the cx-field manifest themselves in observable quantities 
in so far as they affect the fluctuations of the occupation numbers of particles that couple to 
the cr-field, like protons and pions, which interact with a through the effective Lagrangian: 
2G<7 7r + 7r~ +g p a p p. Throughout this paper we will use G = 300 MeV (see Ref. [3]) and g p = 7 
(see, e.g., [9]) as benchmark values. It is important to bear in mind that both these parameters 
and A3 and A4 are all uncertain at the factor of 2 level. 

The contribution of critical fluctuations to the 2-, 3- and 4- particle correlators due to o~- 
exchanges is given by (see Refs. [TU], [3] and [H]) 
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Equations ([T]) - ^ apply to both protons (with g = g p ) and pions (with g = g^ = G/m n ). The 
degeneracy factor d is 2 for both protons and pions since the coupling to the cr-field is both spin 
and charge "blind" . The variance of the fluctuating occupation number distribution is denoted 
by and 7^ is the relativistic gamma-factor of the particle with a given momentum k. 

Mixed pion-proton correlators can also be obtained from ([T]) and Q, with small modifications. 
For example, the 2 pion - 2 proton correlator is given by 

0«7r«p f^ 2 \ \ „2„2t7 U ki ^k 2 V k 3 "k 4 

7T IT P P 

7k a 7£ 2 7k 3 7k 4 



In order to obtain the critical contribution to the cumulants of the particle multiplicity 
distributions using the correlators given above, each momentum index needs to be integrated 
over using V f ci 3 k/(27r) 3 . A cumulant with i protons and j pions will be denoted by Ki P j n . In 
order to remove the V dependence of the cumulants we choose to normalize them as 
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We can see from expressions ([T]) - ([2]) that higher cumulants are proportional to higher powers 
of £ and thus increase by a larger factor near the critical point where £ becomes large. 

In the case of free particles in the classical Boltzmann regime, with no critical fluctuations, 
the fluctuations of any particle number obey Poisson statistics. The Poisson contribution to 
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Figure 1. The ^-dependence of selected normalized cumulants. We only include the Poisson 
and critical contributions to the cumulants. In the left panel we show uj^ p and vary A. In the 
right panel we show other examples of w's with A = 100 MeV. See text for details. 



uji v and Wot is 1, and zero for the mixed cumulants ooipj^. In reality, this 1 of Poisson statistics 
gets few percent contributions from Bose-Einstein statistics, from system size fluctuations, from 
initial state correlations that are incompletely washed out, and from interactions other than 
those with the fluctuations that are enhanced near the critical point. We are ignoring all of 
these non-critical corrections to the Poissonian 1. Existing data on the ratio of the fourth to 
second net proton cumulants, K-A{p-p)/ K 2(p-p)i & t \/s = 19.6, 62.4 and 200 GeV [llj confirm that 
the non-critical corrections to the Poissonian 1 are indeed small. 



2.2. Possible experimental outcomes 

We now illustrate possible experimental outcomes of measurements of the cumulants, assuming 
that the matter produced at the freezeout point of the fireball evolution for some collision energy 
yfs is near the critical point. To start, let us assume that the critical point occurs at /jn = 400 
MeV and that £ max — 2 fm. The simplest ansatz for £(//b) that we have found that incorporates 
the above physics and has roughly the correct scaling behavior from universality of critical 
phenomena is 

f (Mb) = ~ , ^ max ,„ nl/3 , (5) 
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with W(hb) chosen in order to reflect the asymmetry in when approaching the critical 

point from the left (smaller fig) or from the right. We shall define the width A as the width 
in fig within which £ > 1 fm. We also use our benchmark values of the four nonuniversal 
parameters that determine the Ui P j n for a given £, namely g w = G/m n = 2.1, g p = 7, A3 = 4 and 
A4 = 12 and we have allowed for the fact that the chemical freeze-out temperature T decreases 
somewhat with increasing /^b by using an empirical parametrization of heavy ion collision data 
from Ref. [12]. 

Numerical results for some examples of normalized cumulants are shown in Fig. [TJ We can 
clearly see the peak in all the normalized cumulants near the critical point. The results indicate 
that the more protons are involved in the observation measure, the easier it is to identify the 
critical contribution. It is readily apparent that the measurement of these observables in heavy 
ion collisions at a series of collision energies is very well suited to ruling out (or discovering) the 
presence of the QCD critical point. 



Table 1. Parameter dependence of the contribution of critical fluctuations to various particle 
multiplicity cumulant ratios. We have subtracted the Poisson contribution from each cumulant 
before taking the ratio. The table shows the power at which the parameters enter in each case. 
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S. 3. Ratios of cumulants 

In order to locate the critical point, experimental results on multiplicity cumulants will need to 
be compared to the theoretical predictions of the critical contribution to these cumulants. But, 
recall that we had to choose benchmark values for four parameters: g p , g n , A3, A4 and also use 
an ansatz for £(/x#). These parameters are not known reliably or accurately enough to permit a 
quantitative prediction for the effect of the critical point on any one of the cumulants we have 
described. In this Section, we suppose that at some y/s there is experimental data showing 
several of the cumulants significantly exceeding their Poisson values. We ask how ratios of 
cumulants can be used to extract information on £ and the values of the four parameters. And, 
we construct ratios of cumulants that are independent of £ and all the parameters, allowing us 
to make robust predictions for the contribution of critical fluctuations to these ratios. In Table 
[I] we present the parameter dependence of various cumulant ratios (with Poisson subtracted 
before taking the ratios), where we defined 2 A3 — A4 = A' 4 . 

The correlation length £ and the four nonuniversal parameters always appear in certain 
combinations in the multiplicity cumulants and it turns out that we can only constrain four 
independent combinations. The first four rows of Table [T] show ratios that can be used to 
constrain one example of four such combinations. The last five entries in the table show 
parameter-independent ratios. They have no ^-dependence, no dependence on the four poorly 
known parameters, and no ^-dependence. We find that these five ratios are all precisely 1. 

Now let us see how we can use these five ratios in order to locate the critical point. Suppose 
that as you change the center of mass energy yfs of the collisions, experimental evidence for 
peaks shown in Fig. [T] begin to emerge. How do you check in a parameter-independent fashion 
whether the behavior seen in experimental data is consistent with the hypothesis that it is due 
to critical fluctuations? You first subtract the Poisson contributions, and then construct the 
last five ratios in Table I. If the fluctuations seen in this hypothetical data are in fact due to 
the proximity of the critical point, all five of these ratios will be equal to 1, with no theoretical 
uncertainties arising from uncertainty in the values of the parameters. This would be strong 
evidence indeed for the discovery of the QCD critical point. 
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